Pseudogaps in an Incoherent Metal 
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How are the properties of a metal changed by strong inelastic scattering? We investigate this 
question within the two-dimensional t-J model using extended dynamical mean field theory and 
a generalized non-crossing approximation. Short-ranged antiferromagnetic fiuctuations lead to a 
strongly incoherent single particle dynamics, large entropy and resistance. Close to the Mott tran- 
sition at low hole doping a pseudogap opens, accompanied by a drop in resistivity and an increase 
in the Hall constant for both lower temperatures and doping levels. The behavior obtained bears 
surprising similarity to properties of the cuprates. 

PACS numbers: 71.30.+h,74.72.-h,71.10.Hf 



Most of our present understanding of the properties of 
metals is based on Landau's Fermi liquid theory: low- 
energy excitations are coherent quasiparticles with the 
quantum numbers of electrons. This concept has been 
proven to be extremely successful even in systems where 
interactions are very strong, e.g. in liquid '^He or in 
heavy Fermion compounds. In a few classes of materi- 
als, however, most notably the cuprate superconductors, 
the usual Fermi liquid picture appears to break down: 
transport is anomalous, pseudo-gaps open, entropy is 
large and various ordering phenomena appear to com- 
pete with each other ^, ^, |^. This has been taken 
to indicate that new low-energy, long wavelength excita- 
tions like spinous and holons or more conventionally spin, 
charge, current or pair fluctuations play a dominant role 
0. However, a convincing theory based on such scenar- 
ios is still missing. 

In this paper, we want to follow a different and less ex- 
plored route, investigating the possibility that incoherent 
and local excitations dominate as it might happen espe- 
cially at higher temperatures T when strong quantum 
and thermal fluctuations driven by competing interac- 
tions decohere the fermionic excitations. Our starting 
point is the two-dimensional t-J model which describes 
on the one hand the physics of a doped Mott insulator 
and on the other hand the physics of an antiferromagnetic 
(AF) superexchange interaction between nearest neigh- 
bor spins. Long range AF order (in 2d possible only 
at T 0) gets destroyed by a few percent hole doping. 
The resulting spin state is characterized by short range 
AF correlations, and highly incoherent excitations, which 
arc difficult to describe in any conventional many-body 
scheme relying on quasi-particle excitations. The inco- 
herent character of excitations in the cuprates is clearly 
seen in the high electrical resistivity, the large relaxation 
rates for spin and charge, and the large entropy. 

Prominent feature of the underdoped cuprates is the 



pseudogap in the single particle Q and particle-hole spec- 
tra. A plausible explanation for it involves the effect of 
finite ranged fluctuating antiferromagnetic or supercon- 
ducting domains leading to a distribution of local spin 
gaps We will show below that there is a different 
source of pseudogaps arising through nearest neighbor 
exchange coupling and retardation effects in the presence 
of strong magnetic fluctuations. 

Our approximation scheme, based upon the extended 
dynamical field theory (EDMFT) [^, see below, neglects 
most of the longer-range non-local aspects of the prob- 
lem but include the strong inelastic scattering of electrons 
from local magnetic fluctuations. By comparing our re- 
sults to experiments on the cuprates and to numerical 
results for the t-J model, we investigate to what extent 
features like the pseudogap, the large entropy or the Hall 
effect can be described by a strongly incoherent metal. 

Model: The t-J model describes electrons in a 
tight-binding model subject to (i) the constraint of 
at most singly occupied lattice sites (effected by pro- 
jected fermion creation and annihilation operators, = 
c^(l — ti-ct)), and (ii) to an AF spin interaction. 
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where Si = ^^^f a' ^ta'^crcr'^i'^' '^^ ^P^"^ Operator at 
lattice site i, r denoting the vector of the Pauli matrices, 
and t and J couple only nearest neighbors. 

As our goal is to describe the high-T incoherent regime, 
we will neglect most spatial correlations, assuming that 
the self-energy of the electrons is local, Sj:(ix') = T,{uj). At 
the same time, we will keep track of all Li;-dependences 
as we consider a situation where inelastic scattering is 
very strong. We therefore use the so-called "dynamical 
mean field theory" (DMFT) 0, |. Taken as a purely lo- 
cal approximation, DMFT neglects the intersite J term, 
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an important source for inelastic scattering. To include 
its effect, we consider the so-called "extended" DMFT 
(EDMFT) proposed in This approximation is proba- 
bly best visualized [|j by selecting a single site, the "im- 
purity" , out of the lattice. As we want to neglect spatial 
correlations, we can treat the surroundings as an effec- 
tive medium providing a fluctuating environment which 
consists both of electrons and bosonic spin fluctuations 
due to the coupling by t and J, respectively. Local cor- 
relation functions can therefore be calculated by solving 
the following quantum impurity model: 

^^imp = Es. E.clcj:^ + FEs.(ct + h.c.) 

+ Eg^S ■ hq+lY.^Sd{h^+h+_^) (2) 

Here is a projected fermion creation operator 
for the impurity orbital, Ud — ^ad^da- and Sd = 
^'Sa-,cr'd'^fa,(j'da'- The (unrestrictcd) fermion operators 
ci create a fermionic bath, the boson operators hi cre- 
ate a bosonic spin bath (local magnetic field) coupling to 
the impurity degrees of freedom. The effective medium, 
characterized by the fermion and boson energies Ej: and 
LUg' has to be determined self-consistently by identifying 
the single particle Greens function and spin susceptibil- 
ity of the impurity model with the local Greens function 
Goo and local susceptibility xoo of the lattice model 
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Here we use that within EDMFT the Green's function 
GjT and the spin susceptibility xg have a simple k depen- 
dence as both the self-energy and the irreducible 
susceptibility Xir{ito) are taken to be independent of k 
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where efc and Jq are the lattice Fourier transforms of tij 
and Jy, respectively. It follows from (^) that only the 
densities of states N{uj) = Yjj:5{ijj — Ek) and D{lS) = 
Yi^[6{u! — ojq) — 6{lo + ujq)] arc needed and Ek and ojq may 
be assumed to be isotropic in momentum space. Formally 
EDMFT is exact in the limit of infinite dimension, d 
00, if both t and J are scaled proportional to 1/Vd 

The solution of the impurity problem (|^) for given 
N{uj) and D{uj) is difficult. Even for the model without 
the spin boson field hq' (the well-known Anderson impu- 
rity model) dynamical properties can only be calculated 
numerically, e.g. using quantum Monte Carlo, the nu- 
merical renormalization group (NRG), or resummations 
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FIG. 1: The two lowest order contributions to the Luttinger- 
Ward functional <1> and corresponding self-energies. Only di- 
agrams with no line-crossings are taken into account (a gen- 
eralization of NCA). The broken (wavy) line denotes pseudo 
fermion (pseudo-boson) Green's function G/o- {Gt), and the 
solid lines represent the conduction electron Green's functions 
Gca, the curly line the correlator Gha of the bosonic bath. 
Also shown are the self-energies, the local electron Green's 
function Goo and the local susceptibility xoo- 



of perturbation theory like the non-crossing approxima- 
tion (NCA) or the conserving T-matrix approximation 
(CTMA) jo). Unfortunately, all of these methods ex- 
cept for the resummation of perturbation theory and the 
quantum Monte Carlo method are not easily generalized 
to include the spin boson field hq: 

We will therefore employ a conserving approximation 
in which infinite classes of Feynman diagrams are re- 
summed. We are aiming at a level of approximation cor- 
responding to NCA for the usual Anderson model. In 
order to effect the projection onto the sector of Hilbert 
space without double occupancy of the local energy 
level we use a pseudo-particle representation, where the 
singly occupied state is created by pseudo-fermion op- 
erators /o^,c^ =15!: whereas the empty orbital is cre- 
ated by a boson operator &+. Since the local level is 
either empty or singly occupied, the operator constraint 
Q — 6+6 -I- 'Safafa = 1 has to be satisfied at all times, 
which can be enforced by adding a term XQ to the Hamil- 
tonian (j^) and taking the limit A — > 00. The projected 
local electron operators in (||) may then be replaced by 
d^ = b^fa, turning the problem into a many-body sys- 
tem of pseudo- fermions and slave bosons 6, interacting 
with the fermions cr and bosons of the bath. 

ka ^ 

It is essential that in any approximation one stays 
within the physical Hilbert space and does not violate 
the constraint. Therefore, we employ a "conserving ap- 
proximation" specified by a generating Luttinger-Ward 
type functional $ from which all self-energies and corre- 
lation functions are obtained. We employ the simplest 
conserving approximation by considering only the lowest 
order diagrams in V and / (see Fig. |l|) using that the 
effective hybridization V and the exchange field / are 
small compared to the bandwidth St. Within our con- 
serving approximation, both the local Green's function 
of the physical electron and the local susceptibility can 



FIG. 2: The local spectral function plotted versus frequency 
for T = 0.06t and J/t = 0.3 for various hole-doping concen- 
trations S. Inset: The local spectral function for four differ- 
ent J/t = 0, 0.1, 0.2 and 0.3 and T = 0.06t for a doping of 
5 = 0.01. The evolution of a pseudogap of width J is clearly 
visible. 



be calculated as a simple convolution of pseudo-particle 
Green's functions without vertex corrections (Fig. ^). 

It is important to note that our approximation scheme 
does not include the vertex corrections needed to describe 
correctly how the effective interactions with the bosonic 
and fermionic bath renormalize each other. It therefore 
cannot be expected to capture correctly the behavior at 
low T especially close to the quantum critical point where 
AF order is destroyed by doping [|o|. We believe, how- 
ever, that in the incoherent high-T regime which is the 
focus of our study it is unlikely that such vertex cor- 
rections change the physics qualitatively. On a Bethe 
lattice, our EDMFT equations surprisingly are identical 
to those of a t-J model with fully random J, as studied 
within a wstematic large M expansion by ParcoUet and 
Georges [|ll|, who did not find any pseudogaps. We be- 
lieve this to be an artifact of their approximation which 
uses a Bose-condensed slave boson (6) therefore missing 
the incoherent part of the spectral function. 

Results: What happens when inelastic scattering is in- 
creased by switching on a finite J? The effect is strongest 
at small doping as shown in the inset of Fig. Spectral 
weight is pushed below the Fermi energy Ep and a well 
pronounced pseudogap of width J opens. At T = Q.OGt 
the pseudogap closes for (5 ^ 10% as shown in Fig. |^. 

It is tempting to compare our results to experiments in 
the pseudogap phase of the cuprates. One should how- 
ever keep in mind that in the cuprates nonlocal effects 
do play an important role, as is evident from the mo- 
mentum dependence of the pseudogap. Furthermore, it 
is important to stress that we do not see a pseudogap 
in the susceptibility. However, it is interesting to ask 
which qualitative features can be understood as a purely 
local effect. For example, our approximation scheme ex- 



FIG. 3: Entropy as a function of doping for J = 0.3t, 
T = O.lt and T = 0.2t cornpared to results from exact di- 
agonalization (dotted lines) |^| and experiments in LSCO [^] 
at T 0.07t. The triangles mark the doping below which a 
pseudogap starts to open in the spectral function. 



plicitly excludes the possibility that the reduction of the 
density of states at Ep is created by the adjustment of 
the electronic wavefunction to some small magnetic or 
superconducting domains. The observation that a pseu- 
dogap can arise in an incoherent metal with purely local 
correlations is one of the main results of this paper. In 
our approach, the pseudogap opens when the renormal- 
ized chemical potential ReE(u; 0) is pushed towards 
the edge and finally out of the lower Hubbard band jl^] 
by strong magnetic fluctuations: this is only possible in 
an incoherent metal when ImE is sufficiently large. 

How does this physics manifest itself in other physical 
quantities? We calculate the entropy as a crude measure 
for the relevance of incoherent excitation from the free 
energy Q: 
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where the impurity contribution in terms of the pseudo- 
particle spectral functions Af i, is given by e"''^'™'' = 



/ duje 



-0U1 



The entropy S 



—dfl/dT as a function of doping for various T is shown in 
Fig. H. First of all, one realizes that it is rather large even 
at the lowest temperature of T/t = 0.1, an indication for 
strong correlations and a rather incoherent state. The 
overall magnitude of S compares surprisingly well with 
both exact diagonalization and experiments in LSCO (see 
Fig. |3|). Furthermore, it follows the general trend that en- 
tropy is reduced both for large doping where the system 
should become more coherent and at low doping where 
magnetic fluctuations quench the In 2 entropy of a mag- 
netically disordered Mott insulator (for J = the entropy 
increases for (5 — > 0). Interestingly, the drop in entropy 
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FIG. 4: T-dependence of Rh for J = 0.3t. For small doping 
and r — > 0, Rh approaches the value l/{eo5) expected for a 
single hole in a t-J model |l4|. Inset: T-dependence of the 
resistivity multiplied by doping S. The linear T behavior for 
high T flattens for S > 0.1 for T of the order of J. For S < 0.1 
the resistivity drops in the regime where a pseudogap opens. 



towards low doping occurs precisely when the pseudogap 
starts to open in the spectral function (note that in the 
experiment both the opening of the pseudogap and the 
drop in entropy occur at higher doping). 

How is transport affected by the pseudogap? Within 
EDMFT, there are no vertex corrections to the conduc- 
tivities a^x and a^y which can therefore be directly cal- 
culated from the spectral functions ||l^, |l^ . In the inset 
of Fig. ^ the T-dependence of the resistivity is shown 
for various dopings. For high T, p depends linearly on 
T, an effect which is not related to the coupling to the 
bosonic environment as it is also seen in DMFT ||]. For 
small doping, the resistivity is proportional to 1/(5: only 
the holes doped into the Mott insulator can transport 
charge. At the scale of J the resistivity saturates, proba- 
bly due to the strong inelastic scattering from spin fluctu- 
ations. Such a behavior is not observed in experiments in 
the cuprates, possibly an indication that non-local effects 
and vertex corrections are important for transport. Note, 
however, that in the regime where the pseudogap forms, 
i.e. for 6 < 0.1, the resistivity actually shows a clear drop 
which is reminiscent of what is seen experimentally [Q. 

In Fig. U the T-dependence of the Hall constant Rh is 
displayed for various dopings. We find a strongly grow- 
ing positive Rh with decreasing temperature for small 
dopings and an almost flat variation for moderate dop- 
ings. In the limit of small doping and low T the universal 
relation Rh = ^ is approached [ p^ , an indication that 
Luttinger's theorem is not applicable in this incoherent 
regime which cannot be described by moderately excited 
Fermi quasiparticles. Note that the rise of Rh towards 
low T seems to happen in the regime where the pseudo- 
gap opens - in underdoped cuprates a strong increase of 
Rh with falling temperature is observed upon entering 
the pseudogap regime [Q (the drop of Rh close to Tc is 
obviously not included in our theory). In the absence of 



the coupling to the bosonic bath, i.e. within DMFT, both 
the pseudogap |^ and such an upturn [T^ are absent. 

In conclusion, we have investigated the properties of 
a highly incoherent metal close to a Mott insulator sub- 
ject to strong magnetic fluctuation. Even purely local 
magnetic fluctuations change the physics qualitatively at 
small doping: they suppress the entropy and induce a 
pseudogap by driving the chemical potential out of the 
lower Hubbard band. This leads to an increase of the Hall 
constant and a drop in the resistivity. These features are 
reminiscent of the behavior seen in the pseudogap phase 
of the cuprates. This might indicate that some of the 
physics in the cuprates could reflect properties of a highly 
incoherent metal with dominating local fluctuations. An 
interesting open question is to what extent properties of 
such an incoherent metal are universal and independent 
of the details of inelastic scattering mechanisms. 
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